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UNIFORM EXPANSIVITY OUTSIDE THE CRITICAL NEIGHBORHOOD 

IN THE QUADRATIC FAMILY 


ALI GOLMAKANI, STEFANO LUZZATTO, AND PAWEL PILARCZYK 

Abstract. We use rigorous numerical techniques to compute a lower bound for the exponent of 
expansivity outside a neighborhood of the critical point for thousands of intervals of parameter 
values in the quadratic family. We compute a possibly small radius of the critical neighborhood, 
and a lower bound for the corresponding expansivity exponent outside this neighborhood, valid 
for all the parameters in each of the intervals. We illustrate and study the distribution of the 
radii and these exponents. The results of our computations are mathematically rigorous. The 
source code of the software and the results of the computations are made publicly available at 
http://www.pawelpilarczyk.com/quadratic/ 
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We report on some rigorous numerical results concerning expansivity exponents for various 
parameter intervals for the dynamics outside a critical neighborhood in the quadratic family 
of maps. For completeness and in order to motivate our study, we start in Section 1 with 
some general facts on the quadratic family. The reader who is already familiar with the basic 
background may wish to skip directly to Section 2. 

1. Background and basic definitions 

Let I be a compact interval and /: I —>• / a map. It is well known that the dynamics generated 
by the iterates of / can range from very simple to extremely complex. Two basic examples are 
the maps on the unit interval I = [0,1] given by f{x) = :c/10, for which the dynamics is very 
“regular” since all initial conditions converge under iteration to the fixed point at the origin, 
and g{x ) = lOx mod 1, which instead has an infinite number of distinct periodic orbits, an 
infinite number of dense orbits, and other hallmarks of what is generally described as “chaotic” 
or “stochastic” dynamics. It turns out (in a way that can be formalized) that these two quite 
distinct kinds of dynamical behavior are essentially due to the fact that in the first example the 
map is contracting , i.e., the derivative is everywhere smaller than 1, and therefore nearby points 
get closer under iteration, while in the second example the map is expanding , i.e., the derivative 
is everywhere greater than 1, and therefore nearby points are moving away under iteration. 

Things get much more complicated in examples which combine some regions of contraction 
with some regions of expansion. One of the best known and most studied such examples is the 
real quadratic family of maps 

f a (x) =a-x 2 . 

The fixed points of f a are given by solutions of the equation a — x 2 = x which gives x = 
— 1/2 ± \J 1 + 4a/2. Thus f a has two fixed points for every a > —1/4. For a > 0 these lie on 
opposite sides of the origin and letting p a = —1/2 — \J\ + 4a/2 denote the negative fixed point, 
and letting I a = \p a , —p a \ denote the compact interval defined by this fixed point, it is then easy 
to see that for every a € (0, 2] the map 
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is well defined, i.e., f a maps the compact interval I a to itself. The derivative of f a is f' a (x) = — 2x 
and therefore it is easy to check that for positive values of a the derivative at the fixed point 
p a satisfies \f' a {pa)\ > 1 whereas /'(0) = 0. Thus the map is neither fully contracting nor 
fully expanding. Nevertheless, it is easy to check that for small values of the parameter a 
the fixed point q a = —1/2 + y/1 + 4a/2, which lies in the interior of the interval I a , satisfies 
\fa(Qa)\ < 1 and is indeed an attractor for every point in the interior of I a ■ As a increases, 
however, the fixed point q a loses its stability, the proportion of the interval I a in which the map 
is expanding increases, and generally the dynamics gets significantly more complicated. Basic 
numerical simulations obtained by iterating a more or less arbitrary initial condition give rise 
to the well known bifurcation diagram , see Figure [U which was first observed in the 1970’s and 



Figure 1. Bifurcation diagram for the quadratic family. 

generally credited to Feigenbaum and which suggests certain regions of parameters for which 
the dynamics is basically still regular, with most initial conditions eventually attracted to some 
periodic orbit, interspersed by regions of parameter space in which the dynamics is very wild and 
chaotic. Rigorous analytical results over the last couple of decades have gone a very long way in 
explaining and clarifying this bifurcation diagram, and revealed some initially unexpected and 
surprising facts. We restrict ourselves here to the parameter interval 

n = [i.4,2] 

since this is the region in which the bifurcation diagram exhibits interesting phenomena. A first 
and remarkable observation which was already made as far back as the 1940’s by Ulam and von 
Neumann [30] is that the parameter 2 can be thought of as “stochastic” in a precise mathe¬ 
matical sense (more formally, the corresponding map f a admits an ergodic invariant probability 
measure which is absolutely continuous with respect to Lebesgue measure). In 1981, Jakobson 
[6] and then in 1985 Benedicks and Carleson [2], showed that the set C 0 of such stochastic 
parameters is actually “large” in the sense that it has positive Lebesgue measure (see also some 
generalizations in [2311211 23 , [ 2(71 [271 12 , 113 ]). Interestingly, however, it is also “small” in the 
sense that it is topologically nowhere dense, a fact that follows from the remarkable result that 
the set fl” C Ll of “regular” parameters (for which almost every initial condition is eventually 
attracted to a periodic orbit) is open and dense in Ll [H [T5] [TB] (see also generalizations in 
mm)- Thus both regular and stochastic parameters have positive probability in Q. It was 
proved in m that in fact the union U Ll C S2 has full measure in Ll and these are therefore 
the only two dynamical phenomena that occur with positive probability, even though there are 
infinitely many parameters for which different, and sometimes extremely bizarre, dynamics does 
occur. 

On a more quantitative level, some attention has been devoted to the natural question of 
how the measure is shared between the regular and stochastic parameters. Interestingly there 
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seems to be no heuristic argument of any kind suggesting either that most parameters should be 
regular or that most parameters should be stochastic. Recent rigorous numerical estimates in 
m confirm previous calculations in [25] and show that for the related family f a (x) = ax( 1 — x) 
and the corresponding relevant parameter interval [2,4] the set of regular parameters occupies at 
least 10% of all parameters. While the calculations have not been carried out for the quadratic 
family in the form that we consider here, it can be expected that they would yield similar results. 
Parameters in cannot be computed directly since they belong to a Cantor set (which makes 
them theoretically un-decidable, see El), but it is possible to obtain explicit bounds for the 
measure of fi + . Jakobson first set out the theoretical framework and posed the question of the 
possibility of obtaining explicit estimates in f?[ 18]. In [14] it was shown that |fl| > lO -5000 and 
improvements on these estimates have been recently announced by Yu-Ru [5] and Shishikura 
[26], but the bounds obtained are still very small. Some work in progress by the authors of the 
current paper aims to extend these estimates significantly. 

One of the key conditions used in all proofs of the positivity of the Lebesgue measure of the set 
of stochastic parameters is the so-called uniform expansivity outside the critical neighborhood 
which will be defined precisely in the next section. This property holds under certain general 
abstract conditions, but for the quantitative results mentioned above, which include numerical 
bounds on the measure of fi + , it is necessary to have very explicit numerical estimates for the 
expansivity and the size of the critical neighborhood outside of which the expansivity holds. 
The main objective of this paper is to present the rigorous results of such computations for a 
relatively large number of small parameter intervals. 

2. Expansion outside the critical neighborhood 

As mentioned above, the dynamics of a map f: I —>• I is very much influenced by whether / 
is expanding or contracting. Even if / is not expanding we can still formulate the notion of / 
being expanding in certain region of the phase space. 

Definition 2.1. For A > 0 ,/:/—>•/ is A -uniformly expanding outside Ac I if there exists 
C > 0 s.t. for every x G I and every n > 1 with f l (x ) ^ A for all i = 0,1,. .., n — 1, we have 

\(f n y(x)\>Ce Xn . 

Notice that / is trivially A-uniformly expanding outside A if |/ / (a;)| > e A > 1 for all x (j A 
(in which case one could also take (7 = 1), and so in particular any expanding map is expanding 
outside any set A. However, in many situations this property holds even if the derivative of 
/ outside A is very small in some regions. Indeed, a celebrated and important result of Mane 
[18] says that if all periodic orbits of / are hyperbolic repelling (i.e., \(f e )'(p)\ > 1 where i is 
the period of the periodic orbit p) then f is uniformly expanding outside any neighborhood A 
of the critical points (i.e., all the points c such that f(c ) = 0), and even more generally, even if 
/ has attracting periodic orbits then it is uniformly expanding outside any region A as long as 
A contains the critical points and the “immediate basins’’ of the attracting periodic orbits. This 
theorem has multiple consequences and corollaries and is arguably one of the most fundamental 
technical results in one-dimensional dynamics, on which many deep and important subsequent 
results are based. 

Mane’s Theorem is an abstract “existence” result, and a natural question, and indeed the 
question we address and study in this paper, is that of an explicit numerical study of the 
size of the region A and the expansivity exponent A in specific examples, in particular in the 
quadratic family, which we use as a case-study. The computational techniques required for a 
rigorous numerical investigation of this problem were developed in a previous paper [3] . Here 
we further refine those techniques and apply them to a large-scale and systematic study of the 
entire relevant parameter space. 
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In order to formulate our results, we introduce some definitions as follows. In the case of 
maps f a belonging to the quadratic family, with a £ f2, it is known that there can always be at 
most one attracting periodic orbit and that its immediate basin is a neighborhood of the critical 
point c = 0, and that all other periodic orbits are always hyperbolic repelling. This allows us to 
make the following definitions. For each a £ Q, let 

5 a := inf{(5 > 0 : f a is uniformly expanding outside A = (—5, (5)}. 

If f a has a periodic attractor then 5 a > 0, and either 5 a or — 6 a are boundary points of the 
immediate basin of attraction of the periodic orbit. In this case we can define 

X a := sup{A : f a is A-uniformly expanding outside A a = (—5 a ,5 a )} 

and, by the Theorem of Mane mentioned above, we always have A a > 0. If f a does not have 
periodic attractors, then 5 a = 0, i.e., f a is uniformly expanding outside every neighborhood of 
the critical point. In this case we define 

A a := limlsuplA : f a is A-uniformly expanding outside (— e, e)}}. 

e —>0 

In general, we obviously have A a > 0 and, surprisingly, for almost every a £ D + we even have 
A a > o m , meaning that f a is in fact uniformly A-expanding outside every neighborhood of 
the critical point for some A independent of the neighborhood. Notice, however, that f a is not 
A-uniformly expanding for any A if A reduces just to the single critical point itself. 

We now extend these definitions to intervals of parameter values. For a closed non-trivial 
interval oj C D, we let 

:= sup{5 a } = inf{5 > 0 : f a is uniformly expanding outside A = (—5,5) for all a £ u}. 

Since the set of regular parameters with attracting periodic orbits is open and dense, every 
non-trivial interval of parameters must contain some such parameter, and therefore we always 
have 5 U > 0. Given a radius 5 of a critical neighborhood A = (—5, 5), we let 

A(a(5) := inf{A > 0 : f a is A-uniformly expanding outside A = (—5,5) for all a £ u}. 

In the remaining part of the paper, we systematically analyze the parameter space II for 
the quadratic family. We subdivide it into a large number of small subintervals, and obtain a 
collection of rigorous explicit upper bounds > 5 W for 5 U for each of the subintervals u>, and 
lower bounds X u < X U (5 U ), which are expected to approximate A aj (5 a; ). Note, however, that 
X u is, in general, neither an upper nor a lower bound for A w (5 a; ); this is due to the fact that 
the overestimate in calculating 5 W implies that a larger critical neighborhood is considered when 
determining X u , which may increase the chance for computing a larger expansivity exponent 
than for S u . Nevertheless, we prove the following: 

For every a£u, the map f a is X,^-uniformly expanding outside (—5^,5^). 

3. Numerical results and discussion 

We now give plots which summarize the numerical results we obtained. Details about the 
numerical algorithms and procedures will be given in Section 0] together with various remarks 
about the computational issues involved. We emphasize that all the results are absolutely 
rigorous. 

We subdivide the parameter interval fl = [1.4,2] into N = 60,000 adjacent subintervals 
{cujj^Q 1 of essentially equal size; see Section 14.11 for the precise statement and the technical 
details of this subdivision. In order to simplify the notation, we omit the index i when describing 
the computations done for each individual interval Wj. 

For each parameter interval u, we (attempt to) compute a possibly tight upper bound 5 W on 
5 U , as described in Section B~3l The computation was successful for 24,200 out of the 60,000 
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parameter intervals yielding Figure [2] which shows a plot of the value 6^ versus the parameter for 
all intervals of parameters for which such computation was successful. Each of the 24,200 data 
points is plotted as a dot surrounded by a small hollow circle. The circles tend to accumulate in 
some areas, and their superposition may occasionally give the impression of thick (black) dots. 

The lack of success of the calculation for a specific parameter interval u may be due to the fact 
that oj belongs to, or just intersects, some region of regular parameters with attracting periodic 
orbits with relatively large immediate basins, e.g., larger than the maximum size of 0.001 that we 
(somewhat arbitrarily) chose for the maximal allowed radius of the critical neighborhood. This 
is clearly the case, for instance, for the smaller parameters a < 1.55 and for the large obvious 
“window” around a ~ 1.77 corresponding to low period attracting periodic orbits which are 
already obvious from the bifurcation diagram in Figure [I] Although choosing a weaker constraint 
on the radius of the critical neighborhood (e.g., 0.1) would considerably increase the number of 
successful calculations, in a longer perspective we are interested in further arguments regarding 
the measure of parameters for which “stochastic” dynamics occurs, and for this purpose it is 
absolutely necessary that the radius is very small. Moreover, expansivity in the phase space 
is much less interesting if it takes place very far from the critical point only. Therefore, by 
specifying the restrictive constraint on the radius of the critical neighborhood, we focus our 
attention on the results that are important from a wider perspective. 



Figure 2. A possibly small radius S u of the symmetric interval A w around the 
origin for which it was possible to prove that f a is A-uniformly expanding outside 
A w with some A > 0 for all a € u. The circles indicate the pairs (u,5J). 

Parameter intervals for which the computation was successful, on the other hand, indicate a 
good deal of expansion, and therefore very probably a large proportion of stochastic parameters, 
especially if the expansion occurs outside some small critical neighborhoods. One perhaps note¬ 
worthy feature of Figure [2] consists of the several “vertical” rows of dots, occurring for example 
around a ~ 1.9. These seem to suggest the existence of a larger interval, formed by the union 
of a number of the smaller intervals we use in the computation, on all of which we have expan¬ 
sion outside the largest critical neighborhood of radius 0.001, but such that some subintervals 
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continue to show expansion even outside much smaller critical neighborhoods down to intervals 
of radius below 0.0002. Other curious patterns are the short diagonal patterns occurring for 
a ~ 1.7 for which we do not, however, have any heuristic explanation. 

For each parameter interval to for which the computation of 5 U was successful, we compute 
a lower bound for the expansion exponent A outside the critical neighborhood of radius 
following the procedure described in Section 14.21 using a non-uniform subdivision of the 
complement of the critical neighborhood to k = 20,000 intervals (see Section 14.41 for a more 
detailed description of the subdivision and a justification of this choice of the value of k). The 
results of this computation are summarized in Figure [3l 



1.4 1.5 1.6 1.7 1.8 1.9 2 

the parameter a 


Figure 3. Values of X u for which it was possible to prove that f a is A^-uniformly 
expanding outside = (—5^,5^) for all a E uj. The circles indicate the pairs 
(u,X w ) for which the computation was successful. 

Perhaps the main thing to notice here is how reasonably well distributed are the computed 
values of A w , without any particular “jumps” or “thresholds.” We also notice that there are 
quite a few parameter intervals with particularly low values of A w , forming what looks almost 
like a horizontal line very close to 0. It seems that these are significant underestimates of the 
true values of A w , as we discuss in more detail in Section Ed Indeed, we believe that most of 
these estimates represented in Figure [ 2 ] are underestimates due to the fact that improving the 
estimates requires a significant and major increase in required computational cost; we discuss 
this issue more in detail in the next section. 

An obvious question is whether there is any correlation between the computed values of A^ 
and those of <5^; a naive conjecture might be that smaller critical neighborhoods imply a smaller 
expansivity exponent. This is, however, not the case, as shown in Figure Uin which the values of 
A w are plotted against the corresponding size 5 U of the critical neighborhood. Indeed, the picture 
is perhaps most remarkable in the total lack of correlation it shows between these two quantities, 
where both large and small values of X u are essentially equally likely to correspond to relatively 
large or small critical neighborhoods. This is in fact consistent with the known qualitative 
theory of the quadratic family which, as mentioned above, says that for almost all stochastic 
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Figure 4. The computed pairs of values for those subintervals of 

for which the computation was successful, indicated by circles. 

parameters the expansivity constant is uniformly bounded from below for any arbitrarily small 
critical neighborhood. 

We conclude this section with a combined plot of all the three values: the parameter range 
and the corresponding values of 5 and A. This is shown in Figure [5j 



Figure 5. Three-dimensional plot of combined data of the two-dimensional 
plots in Figures El El and |4l 
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4. Computational procedures 

The computation of 6 and A is essentially based on the procedures developed in [3], with some 
minor modifications, including in particular the use of a new and much more efficient algorithm 
in one part of the calculation, see discussion in Section 14.21 below. We clarify selected details 
and briefly outline the methods in this section and refer the interested reader to |3j for more 
comprehensive explanations. 

The source code of the specific software crafted for this paper as well as raw results of the 
computations are available at http://www.pawelpilarczyk.com/ciuadratic/. The computations 
can be run at multiple machines (e.g., at a computer cluster) in a convenient way, using the 
flexible dynamic parallelization scheme introduced in [21j . which is built into the software. 

4.1. Calculation of sub-intervals of Fl. We use double-precision floating-point numbers fol¬ 
lowing the ANSI/IEEE 754 standard, further called representable numbers , which are available 
from the level of the C+-1- programming language, and implemented at the hardware level 
in modern processors. For a real number a, let a denote its representable approximation. 
In order to subdivide the parameter interval [1.4,2] into N = 60,000 adjacent subintervals 

Zj 1 of equal size (up to rounding), we first select the subdivision points. Namely, for each 
i G { 0 , ..., N}, define 

9 q , (i/gcd(i,AQ)(2 - L4) 

*' N/ gcd(«, N) 

where gcd(i, N ) denotes the greatest common divisor of i and N, and all the operations in 
the formula are done in double-precision floating-point arithmetic, with rounding the result of 
each operation to the nearest representable number (this is hardware-supported). Then we set 
Ui := The purpose of truncating the fraction i/N by the gcd (i,N) is that we would 

like to make sure that the endpoints of the intervals computed for a finer subdivision, e.g., for 
N = 120,000, agree with those for a coarser one, which need not be the case otherwise, due to 
differences in rounding. 

In order to avoid notation overload, in what follows, whenever explicit values are mentioned for 
endpoints of an interval w, e.g., [1.9999,2], we always mean their representable approximation. 

4.2. Computation of a lower bound for the expansion exponent A. Given an interval 
uj C Ft and a real number <5 > 0, we compute a lower bound for A > 0 such that f a is uniformly 
expanding outside A := (—6, 5) for all a € u>, using a graph approach introduced in [3] with a 
non-uniform subdivision of the complement of the critical neighborhood into a given number k of 
intervals. However, we introduce one major improvement. Namely, we use a new algorithm [22] 
for computing the minimum cycle mean in a directed graph, instead of original Karp’s algorithm 
[9], which was applied in [3]. This tremendously reduces the memory usage, from 0{k 2 ) to 0{k). 
For example, the previous algorithm required some 400MB RAM for k = 5,000, and with the 
new approach the memory usage drops down to some 12MB for the same graph. This difference 
is more profound for large graphs: Over 50GB RAM required for k = 50,000 now reduces to 
some 20MB. This improvement allows us to use much larger values of k than it was possible 
previously in [3], and shifts the bottleneck from memory requirements to time constraints; see 
the discussion in Section sa 

The basic idea of our approach to the computation of a lower bound for A is to reduce the 
problem to one of bounding the mean weights of paths in certain weighted digraphs (directed 
graphs) related to the map / under consideration. We let G = (V, E, w) denote a weighted finite 
digraph , where V denotes the finite set of vertices, E C V x V is the set of edges, and w : E —» M 
is the weight function. A path is a nonempty finite sequence of edges F = (e\,..., e n ) such that 
e j = (vj,Vj) € E and vj = Vj +1 . The path T is called a cycle if v\ = v\. The weight and mean 
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weight of a path T = (ei,..., e n ) are defined by 

W(T) = ^ j w{e j ) and W(r) = 

3 =1 71 

respectively. 

In the setting of a quadratic map f a as above, and a critical neighborhood A, we consider a 
collection of intervals X = {Ij \ j = 1 with pairwise disjoint interiors and which cover 

I\ A, and say that the weighted digraph G = ( V , E, w) is a representation of f on I\ A provided 
that: 

(a) V = X U {cl A}; 

(b) {e = (Ii,/ 2 )€lx7|/(/ 1 )n/2^0}c£ 

(c) For each e = (h,h) £ E, w(e) < inf { log |X>/(x)| : x £ I\ n 

Observe the following straightforward relationship between the weight of a path and the 
derivative along points whose orbit is described by the path. Given a point x £ I \ A and a 
path T = (ei,..., e n ) such that ej = {Ij-i- Ij) and fi(x) £ Ij for all j = (),..., n, we have 

n— 1 

(1) log|£>/"(*)| =^log\Df(fi(x))\ >W(T). 

3=0 

The representation of the map and its derivatives to a weighted graph, and in particular the 
bound (HD, reduce the problem of determining expansion estimates to the computation of mean 
weights of certain paths. Specifically, the minimum mean weight of any cycle in G provides a 
lower bound for the exponent A of interest. This is the quantity computed by Karp’s algorithm [9] 
or its improvement [22]. 

Some remarks were made in [3] concerning the effect of using partitions with larger or smaller 
numbers of elements, or partitions with elements of variable size. As a result of those discussions, 
and a few additional tests, we choose a non-uniform partitioning of the interval I u := [j a£uj Ia, 
as an apparently more effective one. 

Note that the computations are done using rigorous numerics (with controlled rounding di¬ 
rections), and are conducted for the entire range u of the parameters a at a time (by taking 
lower or upper bounds on the corresponding values, where appropriate). 

4.3. Computation of an upper bound on 5 U . Given an interval ui, a possibly tight upper 
bound for is computed in the following way. For 6 > 0, let Af, denote the lower bound for 
the expansion exponent outside (—5,5), computed with the procedure described in Section fT2l 
for the parameter interval u, using a coarse partition of k = 1,000 subintervals. Starting with 
an apriori chosen 5o := 0.001, we first compute Af, 0 . If Af, 0 < 0 then the computations for oj are 
considered to fail, and no suitable upper bound on 5 U is reported. Otherwise, a possibly small 
5i £ (0,5o] is found for which Af, 1 > 0, using 20 steps of the bisection method applied to [0, 5o]. 
This number is further called 5 U and constitutes a rigorous upper bound for 5 U . 

At this point we remark that in spite of the successful computation of a positive expansivity 
exponent for 5 W using the coarse partition of k = 1,000 subintervals, in some rare cases it may 
happen that A w computed using a finer partition turns out to be negative. This may be due 
to differences in the subdivision of the interval, because of using a non-uniform partition, and 
should be treated as a numerical artifact. If this happens, we also treat this situation as a failure 
in the computation of 5 U and A w . 

4.4. The choice of the partition size k. As it was already discussed and shown in [3], the 
larger the number of partition elements k > 0 is chosen, the higher the lower bound for the 
expansion coefficient A is found. However, the time complexity of the algorithm for computing 
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this bound is 0(k 3 ), which practically means that the cost increases considerably with the 
increase in k. 

In order to find a reasonable number k for our comprehensive computations, we conducted a 
test with u = [1.9999, 2] and a selection of different values k, ranging from 1,000 to 40,000. The 
computed values of A w as a function of k (note the logarithmic scale at the k axis), as well as 
the computation time (note the logarithmic scale at both axes) are indicated in Figure El 
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Figure 6. The exponent A w (the plot at the left-hand side) and the time tk of 
computation of the exponent X u (the plot at the right-hand side) as a function 
of partition size k, computed for oj = [1.9999,2], 
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Figure 7. The smallest value of k E {1,000,2,000,5,000,10,000,20,000} for 
which the computed X u is above 0.001. Circled dots at the top edge of the figure 
indicate the cases in which the value of X u is below this threshold for all the values 
of k in the selection. The number of dots in each row is indicated in parentheses 
at the right-hand-side of the diagram. 
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As one can see, the computed exponent A w is essentially 0 up to a certain threshold, and 
remains nearly constant afterwards. It seems that taking too low a value for k does not allow 
to see the expansivity, which is then recovered after a certain threshold. We conducted the 
computation of A w for N = 60,000 and a few values of k between 1,000 and 20,000; see Figured 
This experiment shows that the discussed threshold is different for different parameter intervals. 

Out of the 24,200 subintervals for which the computation of the expansivity exponent was 
successful, a value A w > 0.001 is encountered in all but 335 cases when k reaches 20,000, which 
constitutes over 98.6% of all the cases. Taking this into consideration, and also bearing in mind 
the cost of the computations (which increases nonlinearly with the increase in k), we decided 
to use k = 20,000 for all the plots discussed in Section [3j Indeed, this value of k seems to be 
beyond the threshold that allows the computation of the seemingly correct exponent A w for a 
vast majority of intervals oj, yet the computation cost is reasonable and takes about 4 days on 
a relatively modern server capable of running 32 processes at its full speed. Therefore, it seems 
that there is no point in increasing k any further. This justifies our choice of k. 
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